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RR LYRAE - THEORY vs OBSERVATION 
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ABSTRACT 

The luminosities, effective temperatures and metallicities that are derived empirically by Kovacs and 
Jurcsik from the light curves of a large number of globular cluster and field RRab and RRc stars are 
compared to theoretical RR Lyrae models. The strong luminosity dependence of the empirical blue and 
red edges (LogL vs. LogT c g- diagram) is in disagreement with that of both radiative and convective 
models. A reexamination of the theoretical uncertainties in the modelling leads us to conclude that the 
disagreement appears irreconcilable. 

Subject headings: stars: variable, stars: oscillations, stars: RR Lyrae, stars: horizontal-branch, stars; 

abundances, stars: distances, stars: atmospheres, stars: fundamental parameters, 
globular cluster 
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1. INTRODUCTION 

The recent work of Kovacs & Jurcsik (Kovacs & Jurcsik 
1996 [KJ96], 1997 [KJ97], Jurcsik 1998 [J98]) proposes an 
almost purely empirical method of extracting the absolute 
magnitudes, colors and effective temperatures (T c ff) and 
metallicities directly from the observed periods and light 
curves of RR Lyrae stars. The only theoretical input ap- 
pears in the transformation from My and V - K to L and 
T e ff via Kurucz's static model atmospheres. The very em- 
pirical nature and potential usefulness of the approach has 
attracted a great deal of attention from observers. 

The most recent work (J98) compiles and analyzes a 
sizeable set of observational data of RR Lyrae cluster vari- 
ables and RR Lyrae field stars. This study thus includes 
RR Lyrae stars with metallicities ranging from Z—0. 00001 
to almost solar Z=0.020. The end products of her analysis 
that we are most concerned about here are Log L-Log X^f 
and LogL-LogP plots, and the byproducts which are re- 
lations between luminosity L, mass M and metallicity Z. 

The J98 Log L-Log T e g and Log L-Log P data are shown 
in Figure 1 as small dots, circles for RRab and triangles 
for RRc stars. The figures indicate well defined slopes for 
the fundamental and overtone blue edges and red edges, 
all four of which have approximately the same values, al- 
though there may be a slight broadening of the instabil- 
ity strip with luminosity. Hereafter we call this slope 2 
(= A Log Lj A Log T cff ). 

J98 pointed out that no model calculations can explain 
the strong dependence of the temperature on the lumi- 
nosity. In §2 we first confirm that indeed the slope 2 of 
this empirical LogL vs. LogT c ff diagram is in irreconcil- 
able disagreement with that of radiative models. Next, in 
§3 we show that the inclusion of turbulent convection in 
the models only shifts the LogL vs. Logoff line, but does 
not change the slope. The discrepancy is therefore very 
basic because it is already inherent in models that are as 
fundamental as purely radiative ones. In §4 we reexamine 
the uncertainties of both radiative and turbulent model 



calculations and conclude that the theoretical 2 is very 
robust and cannot be changed very much without intro- 
ducing new physics. In §5 we discuss the theory that goes 
into the empirical relations of Jurcsik, in particular the in- 
fluence of the use of the static Kurucz color - temperature 
transformation. We conclude in §6. 

2. RADIATIVE MODELS OF RR LYRAE STARS 

In the following we examine how well RR Lyrae mod- 
els agree with the J98 data. We first examine radiative 
models, and thus limit ourselves to the vicinity of the blue 
edges. 

Figure 1 shows a Log L-Log T c ff plot on the left. On 
these we superpose the linear fundamental and first over- 
tone blue edges of radiative RR Lyrae models. We recall 
here that among the linear edges of the instability strip 
only the overtone linear blue edge and the fundamental 
linear red edge coincide with the observable edges of the 
instability strip. Nonlinear dynamical effects shift the ob- 
servable fundamental blue edge and the overtone red edge 
to the red and to the blue, respectively, compared to their 
corresponding linear edges, by up to several 100 K (e.g., 
cf. Buchler 2000). Strictly speaking, for radiative models, 
it is only the linear overtone blue edge that is relevant, 
but the linear fundamental blue edge can give a rough in- 
dication of the slope 2. Furthermore, the red edges are 
determined by convection, and radiative red edges are not 
relevant and are therefore not shown in Fig. 1. 

The solid lines have been computed for models with 
M=0.65, X=0.75 and Z=0.001. The dotted lines repre- 
sent models with M and Z from Jurcsik's {L, M, Z} re- 
lations, the dashed lines those with M from Jurcsik, but 
with Z=0.001. Finally, the large open circles are mod- 
els from Tuggle and Iben (1972) for M=0.6 and Z=0.001 
(with the old Los Alamos opacities). They have essentially 
the same 2 as ours which are computed with the OPAL 
opacities (Iglesias & Rogers 1996) merged with the low 
temperature Alexander-Ferguson (1994) opacities. 
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Fig. 1 : Left: LogL-LogT c g plot: RRab (filled circles) and RRc (triangles) stars from Jurcsik ([J98] and priv. comm.); Theoretical 
fundamental and first overtone blue edges for radiative models: solid lines: fixed mass and Z; dotted lines have M and Z from 
Jurcsik's relation; dashed lines have M from Jurcsik with fixed Z=0.001; Tuggle & Iben's overtone and fundamental blue edges are 
reported as open circles and squares, respectively. Right: Same data in a LogL-LogP plot. 



Again, because we are computing radiative models we 
do not expect the location of the blue edge to be in per- 
fect agreement with the edge of the instability strip, but 
we note a very large discrepancy in the slope S. 

Figure 1 also displays a LogL - LogP plot on the right. 
Only the slope of the constant mass models is in almost 
acceptable agreement with the J98 data, but those calcu- 
lated with the J98 {M, L, Z} are in strong disagreement. 

In Figures 2 we display again a LogL - LogT e jf plot 
which shows the effect of composition on the location of 
the linear blue edges. The (radiative) models were com- 
puted with M=0.65. Here the solid lines have X=0.75, 
Z=0.001, the dotted lines have X=0.75, Z=0.004 and the 
dashed lines have X=0.70, Z=0.001. 

The location of the blue edge displays very little sensi- 
tivity to either mctallicity Z or helium content Y, within a 
reasonable range of values. Adjusting the helium content 
or mctallicity does not provide a resolution of the slope 
discrepancy. We describe further tests in §4. 

We conclude this section by noting that purely radiative 
RR Lyrae models are in severe disagreement with Jurcsik's 
Logi-LogT e ff data. 

3. CONVECTIVE MODELS OF RR LYRAE STARS 

The addition of time-dependent turbulent convection in 
the models has led to some spectacular successes compared 
to radiative models, in particular, in finally predicting dou- 
ble mode behavior both in RR Lyrae (Feuchtinger 1998) 
and in Cepheids (Kollath, Beaulieu, Buchler & Yecko 
1998). One would hope therefore that the inclusion of 
turbulent convection might also remove the discrepancy 
between observations and theory described in the previ- 
ous section. 



Yet, in that respect all of our modelling efforts with tur- 
bulent convection have proved in vain, despite the flexibil- 
ity afforded by the 8 free order unity (a) parameters that 
the turbulent convective equations contain (e.g., Yecko, 
Kollath & Buchler 1997, Buchler& Kollath 2000). We 
find that turbulent convection can shift the blue edges, 
but cannot produce the differential effect with respect to 
luminosity that is required to give the right slope S in the 
LogL vs. Log T c ff plots. In Figure 3 we show the results 
for three different combinations of these parameters. 

We note that different codes and slightly different 
recipes for convection give essentially the same theoreti- 
cal slopes. For example, even though no linear models are 
computed, the nonlinear hydrodynamical models of Bono 
et al. (1997) indicate a similarly steeper slope than that of 
J98. 

4. DISCUSSION OF THE MODELING ASSUMPTIONS 

In this section we discuss the effects of several of the ap- 
proximations that are inherent in the numerical modeling 
to see how robust the theoretical Log L-Log T ff slope is. 

4.1. Equation of State 

Our code uses as equation of state a simple iteration of 
the Saha equations (e.g. Stellingwerf 1982). This equation 
of state is very similar to the OPAL (Rogers et al. 1996) 
and the MHD (Dappen et al. 1988) equations of state and 
we do not believe that the tiny differences can be respon- 
sible for the discrepancy between the models and the J98 
data. 
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4.2. Radiative Transport 

Our radiative code uses a standard equilibrium diffusion 
approximation for the radiative transport, i.e., 



L=(47r i^f A aT 4 (1) 
■in am 

in which one uses for the opacity k the Rosseland mean 
(e.g., Mihalas & Mihalas 1984). Eq. | has several short- 
comings in the region above the photosphere. First, it 
implies an Eddington factor of Je = 1/3, which is not cor- 
rect in the optically thin outer region where Je approaches 
0.4-0.5 (Feuchtinger & Dorfi 1994, Fig. 3). Second, in this 
regime, the opacities should be higher than those given 
by the Rosseland mean (Alexander & Ferguson 1994). Fi- 
nally, for the computation of the periods and growth rates 
it would be more appropriate to linearize the radiation 
hydrodynamics equations rather than the equilibrium dif- 
fusion equation. The magnitude of the errors introduced 
by the first effect can easily be estimated. 

One can approximate the effect of the Eddington fac- 
tor in Eq. |l|. by replacing 3k by k//e- The Eddington 
factor also appears in the radiation pressure which be- 
comes p ra d — > 3/ePrad- We can disregard sphericity ef- 
fects in the momentum equation contained in the term 

(Prad - fE Grad)/?". 

To simulate the effect of an increasing we assume that 
for T < T e s, Je increases smoothly from 1/3 to 1/2 with 
decreasing temperature. The results confirm our hunch 
that an increase of /e in the very outer region has very 
little influence on the growth rates and almost none on the 
periods. More importantly it has no differential sensitivity 
to luminosity. 
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Fig. 2 : LogL-LogT c ff plot: RRab (filled circles) and RRc (tri- 
angles) stars of Jurcsik; Radiative blue edges for M=0.65; solid 
line: X=0.75 Z=0.001, dotted line: X=0.75 Z=0.004, dashed 
line: X=0.70 Z=0.001. 
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Fig. 3 : LogL-Logr c ff plot: RRab (filled circles) and RRc (tri- 
angles) stars: Convective models: Solid, dotted and dashed 
lines are for three different combinations of the a parameters 
in the convective model equations, chosen to give a reasonable 
width for the instability strip. 

We have artificially, and somewhat arbitrarily, increased 
the opacity from the vicinity of T e $ outward. We have not 
found an appreciable differential effect with T c ff, and have 
therefore not pursued this avenue with more detailed mod- 
elling. 

Finally, the insensitivity of the slope to Je makes us 
believe that the linearization of the radiation hydrody- 
namics equations instead of the common equilibrium dif- 
fusion equation (requiring a serious coding effort) would 
also make little difference. 

We conclude that the discrepancy between the models 
and the observations is not due to an inadequate treatment 
of the radiation transport. 

4.3. Composition 

Could compositional make-up be more important than 
expected? We have computed a number of RR Lyrae mod- 
els with various combinations of Y and Z, but find that the 
effect on S is negligible. Next we have artificially increased 
the abundances of the easily ionizable light elements such 
as Mg and Na in the OPAL opacities. Again this has such 
a negligible effect on 3 that the inclusion of this data would 
unnecessarily clutter Figure 2. 

An inhomogeneous composition is unlikely to exist in 
RR Lyrae envelopes because of convection, but even if ex- 
isted our tests with various changes in composition make 
us doubt that it would resolve the discrepancy. 

4.4. Rotation 

Could rotation be responsible for the discrepancy? In 
order to estimate the magnitude of the effect of rotation we 
have included a spherical pseudo-centrifugal acceleration 
u> 2 r in the equilibrium model and in the computation of 
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Fig. 4 : Left: T e g-(V-I) plot and right: (V-I)-Mv plot: Theoretical color variations for two nonlinear full amplitude RRab models 
with different luminosities of L = 40 Lq (dashed line) and L = 52 Lq^uII line). Crosses indicate the time steps of the nonlinear 
calculations. See text for details. 



the periods and growth rates. We find that a rather short 
rotation period, of order of a few days, would be neces- 
sary to have an impact. Furthermore the disagreement is 
worsened by rotation because the slope steepens. 

We conclude that moderate rotation rates cannot be the 
cause of the discrepancy. 

4.5. Evolutionary Effects 

The blueward moving evolutionary paths turn around 
at some T e j and then move upward (e.g., Dorman 1992, 
Lee, Demarque & Zinn 1990). If this happened inside the 
linear instability strip, then the leftmost blue edge, that 
of the overtone mode, would be defined by the topology of 
these paths, rather than by the pulsational stability of the 
models. However, the fact that the Jurcsik blue and red 
edges of the fundamental and the overtone are all essen- 
tially parallel eliminates this possibility as an explanation 
for the shallow observed slope. Moreover, for some of the 
clusters, full photometry of the horizontal branch exists, 
and there is no indication of any gap next to the RR Lyrae 
region. Clearly there are stars next to the empirical red 
edge for which the amplitude (if not zero) is less than the 
observational limit. 

4.6. Conclusion 

The theoretical S slope is very robust with respect to the 
uncertainties or approximations inherent in our code, but 
is in disagreement with the Jurcsik empirical relations. We 
have therefore been led to reexamine the latter for their 
robustness, especially since they also make use of some 
theory, viz. color to temperature transformation. 



5. DISCUSSION OF THE EMPIRICAL RELATIONS 

Because of the strong disagreement between Jurcsik's 
processed observational data and the predictions from 
something as basic as radiative models, it is worth while 
to reexamine some of the most uncertain points in the re- 
duction of the observational data. 

The observational material consists of RRab periods and 
lightcurves. KJ96, JK97 and J98 deduce the color, the vi- 
sual magnitude M v and the metallicity Z ([Fe/H] in the 
form of linear fully empirical relations 



M v \ 
V-K 
[Fe/H] J 
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I P \ 
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hi 

V04l7 



(2) 



The connection with LogT e jf and with LogL was then 
made through static envelope models of Kurucz (1993, cf. 
J98 for details). This is the only theoretical input into the 
otherwise empirical relations. 
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5.1. Color - T c g Transformation 

J98 uses the static envelope calculations of Kurucz's AT- 
LAS code (for details cf. J98), to make the transforma- 
tion from average color to average T e g. The T e g however 
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Fig. 5 : Left: RRab Period - Radius plot; open circles: Baade-Wesselink radii, and crosses: same stars from J98; The solid and 
dashed lines are linear regressions through the Baade-Wesselink and the Jurcsik data, respectively. Right: thick lines: radiative 
models with Jurcsik L, M, Z relation; thin hues: radiative models with constant mass and luminosity; filled squares: same stars 
from J98; 



varies considerably during one pulsation cycle, compared 
to the width of the instability strip, and the question arises 
whether the correct mean color and mean temperature are 
still related through the static Kurucz relations. To check 
the validity of this assumption we have computed the spec- 
tra of two full amplitude RR Lyrae models. 

We investigate the theoretical color variation for two RR 
Lyrae models with different luminosities (L = 40 Lq and 
L = 52 L Q , both for M = 0.65 M©, T cS = 6500 K, X = 
0.76 and Z=0.001). Nonlinear full amplitude models are 
taken from Feuchtinger (1999) and the color variations in 
V and I are computed as outlined in Dorfi & F euch tinger 
(1999). The results can be inferred from Fig. 4.2 which 
displays the color difference V - I as a function of the 
effective temperature (left panel) and the corresponding 
color-magnitude diagram (right panel), both over one os- 
cillation cycle. Full and dotted lines refer to the L = 40 
Lq and L — 52 L Q models, respectively. The differential 
effect with L is seen to be small. Dynamical atmospheric 
effects thus cannot account for the discrepancy. 

There remains the question of the consistency of the 
static Kurucz atmospheres with our frequency dependent 
calculation. For that purpose we have checked the color - 
temperature transformation against the Kurucz tables at 
selected points on the pulsation cycle. The agreement is 
quite satisfactory, and any discrepancies are irrelevant for 
our purpose here. 



5.2. RR Lyrae Period - Radius Relation 

The radii of RR Lyrae stars can be obtained indepen- 
dently either from a Baade-Wesselink approach or from 
the J98 empirical L and T G ff- In Fig- on the left, we 
have plotted as open circles the Baade-Wesselink radii of 
the 15 RRab stars of Jones et al. (1992) (we have omitted 
RR Leo) . Superposed as crosses are the radii for the same 
stars derived from the J98 L and T c g. The agreement is 
remarkably good, but worsens for the longer periods. The 
solid and dashed lines represent linear regression lines for 
the BW and for the J98 period - radius relations. 

We have also derived the radii of the additional J98 stars 
that were not in the Jones sample and show them as small 
dots in the righthand figure. The thin lines represent con- 
stant mass and constant luminosity radiative models, and 
the thick lines radiative models that obey the J98 L, M, 
Z relations Since the radiative models do not provide a 
red edge, the upward extent of the thin or thick lines is 
not significant. The latter models are seen to give better 
agreement with the swarm of the J98 empirical radii. 

We conclude that the agreement between the empirical 
J98 and Jones' Baade-Wesselink radii is remarkably good, 
and so is the agreement that can be achieved with theo- 
retical models. A consideration of the period - luminosity 
relations therefore does not give us any clues as to the 
origin of the slope discrepancy. 
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6. CONCLUSIONS 

We have shown that the Jurcsik relations lead to RR 
Lyrae radii that are in good agreement with both Baade- 
Wesselink radii and with theoretical radii. On the other 
hand, there exists a strong discrepancy between the slope 
of the theoretical Log L- Log T c g relation and the slope of 
the empirical Jurcsik relation. 

We have reviewed the physical and numerical uncer- 
tainties that enter the theoretical calculations. The dis- 
crepancy exists already at the level of purely radiative 
modelling. We have further shown that our model equa- 
tions for turbulent convection do not alter the slope of 
the Log L-Log T c g- relation despite the large number (8) of 
adjustable parameters (a's) that we have at our disposal. 
Consequently, apart from the uncertainties inherent in a 
ID recipe for turbulent convection, it appears that the dis- 
crepancy is not caused by a deficiency of the theoretical 
models. We note that the same slope discrepancy is im- 
plicit in older and more recent independent calculations 
(Tuggle & Ibcn 1972, Bono et al. 1997). 



The derivation of the empirical relations makes use of 
a static Kurucz color - T c $ transformation. By computing 
the behavior of color versus temperature over the pulsa- 
tion cycle we have shown that the use of static envelopes 
is in fact a very good approximation, and therefore cannot 
be the culprit. 

Finally, the shape of the evolutionary tracks, through 
their potential avoidance of certain regions of the instabil- 
ity strip cannot be responsible either. 

At this time the origin of this disturbing discrepancy 
constitutes an unsolved puzzle. It remains to be seen 
whether an improved treatment of turbulent convection 
or more complete observations will resolve the difficulty. 
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